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CN ■ Abstract 

o 

■ Results of numerical simulations of fusion rate d(d,p)t, for low-energy deuteron 

i beam, colliding with deuterated metallic matrix (Raiola [1,2]) confirm analytical 

O I estimates given in Ref. [3] (M. Coraddu et al, this issue), taking into account 

' quantum tails in the momentum distribution function of target particles, and predict 
an enhanced astrophysical factor in the 1 keV region in qualitative agreement with 

Q I experiments. 

> 

X 

' 1 Introduction 



Significant divergence from theoretical predictions of non-resonant fusion cross- 
section at low energies of incident deuteron particles has recently been ob- 
served in experimental works [1,2]. At energies of the charged deuteron beam 
less than 5 keV, colliding with deuterated metallic matrix, a great enhance- 
ment of the fusion cross section takes place, compared to theoretical evalu- 
ations. The hypothesis of ions interacting via screened potential inside the 
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metallic target was proposed in [1], but an unrealistically large screening po- 
tential was necessary to obtain reasonable agreement with the observed data. 
It is obvious that problems arising from experimental results need refining the 
theoretical models actually used. 

As it was shown in papers [4,5], the role of quantum corrections to the par- 
ticle momentum distribution is quite important. A significant deviation from 
Maxwellian distribution function appears. At large values of momentum, due 
to quantum corrections, the tail of the distribution has a power-law asymptotic 
behavior instead of exponential, resulting in modified reaction rates. In partic- 
ular, this effect leads to non-exponential temperature dependence of inelastic 
process rates at relatively low temperatures and high densities. 

In this paper we present results of numerical calculations of reaction rates for 
conditions of the interacting particles in the beam and in the target such as 
those of the experiments [1,2]. 



2 Numerical modeling of reaction rates for experimental condi- 
tions 

The state of the system defined by the generalized distribution function over 
energy and momentum may be presented in the factorized form: 

F{E, s) = n{E)a{E — s), where s = /2m is the particle kinetic energy. The 
reaction-rate constants of the inelastic process between two particles, named 
"a" and "b", may be presented in a more general form by the integral (see 
Ref. [6], A. Starostin et al.,this issue) 

oo oo 
NaN,Ka,^ST^C J dEaJdpa I dE.JdpJdqX 
— oo — oo 

n{Ea){l-n{Ea + Qa-uj))aa{Ea-ea)x 
n{Eb){l - n{Eb + oo + Qb))ab{Eb - eb)x 

a'^{Ea + Qa-uj- ep^-^)a'b{Eb + uj + Qb- ep^+^)a{ep) 

where Ea , Pa are the energy and momentum of the "a" particle, Sp = p^ /2M 
is the relative kinetic energy in the center of mass, M is the reduced mass of 
colliding particles, C is a normalization constant, defined from comparison of 
the expression calculated by (1) and known results at high temperature and 
low density. We have: 

rUbPa - rriaPb 

P= , ■ 

rria + nib 
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The expression for the population number n{E) depends on the statistical 
distribution of the system. For the purpose of this work we must consider that 
deuterons, which are bosons, have the distribution 

n{E)^{eM{E-pi)/T]-l)-' 

where fi is chemical potential. For non-ideal plasma, the spectral dependence 
of the distribution function, defined by Lorentzian profile, is: 



a{E,p) 



TT 



{E-e,-A{E,p)Y + jf{E,p) 



(2) 



In Eq. (2) 7 is the width, A stands for the energy shift due to atom-matter 
interaction. The line width is given by 7^ = NhaaVa, where cTq = Tre^/s^, 
Va = y^2£a/ma, N is the concentration of scattering centers. For ideal plasma 
conditions, i.e. when the density decreases, the width '~f{E,ep) , the 
function a{E,p) becomes a delta-function. The cross-section dependence on 
kinetic energy may be given in the form: 



where ri{ep) is the Sommerfeld parameter. The astrophysical factor S is weakly 
dependent on the kinetic energy. 

The influence of the screening potential Ue on the reaction rate, due to the 
effect of metallic electrons, may be taken into account by adding Ue to the 
collision energy: 

<7{Ep) = (^oi^p + Ue) . (4) 



Numerical calculations of the reaction rates were performed in accordance 
with the above model, for conditions close to the experimental ones: target 
particles concentration Na — 5 ■ 10^^ cm~^, and interacting particle masses: 
TTia = rrib = 2 amu. Fusion reactions were considered between target particles 
of kind "a" and beam particles of kind "b" . In the expression of the width 
we used as concentration the concentration of the scattering ions in the 
metallic matrix. 

Taking into account high dimension, the computation of the integral (1) can 
be performed using Monte Carlo method. The kinetic energy distribution of 
target particles was taken at temperature T = 2.44 • 10~^eV, while the "b" 
particles were taken with the beam energies. It is interesting to note, that for 
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the case of a beam of mono-energetic particles the reaction rate computation 
can be reduced from expression (1) to the more simple integral: 



NaK' j dEa j dpna{Ea)a{Ea - Ep, ep)J -^a{ep)dep . (5) 





Within the framework of such model, for the case of ideal plasma, when we can 
neglect the wings of the Lorentzian profile in the integral (5), the expression 
for the reaction rate can be further simplified: 



N^K^ = nX j dE^Jdp J na(E,)5(E, - e^)J ^a{ep)dep ~ 

° oo (6) 
J na{ep)epa{ep)d£p . 



oo 




The influence of the distribution wings on the reaction rate value can be 
obtained by comparison of the computation results of the two expressions 
(1) and (5). We can also compare such results with the calculated reaction 
rate Ki — aV, using expression (3). After such comparison it is possible to 
estimate the astrophysical factor S and the deviation of theoretical predictions 
from experimental data. 

To estimate the influence of momentum distribution tails on the reaction rate 
and the difference with the Maxwellian case, it is necessary to take into account 
the finite width of the Lorentzian profile (2). As was shown in [5] and in [6], 
the main result of quantum corrections is that the momentum distribution 
function has asymptotically a power-law tail: 

oo 

fisp) = C' J dE,a{E, - Ep, Ep) - exp{-Ep/T) + Ca{T)/El . (7) 



Using such decomposition in equation (5), it is possible to calculate the reac- 
tion rate, taking into account non-Maxwellian distribution function: 

oo 

K,^C,jdEj{Ea)^a{Ep) . (8) 


The results of such calculations are shown in the table 1. 

The reaction rate constants, calculated by different models, agree among them- 
selves at beam energies above 2 keV. Decreasing energy in the range between 
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Table 1 

Comparison of the reaction rates {ov) as function of the energy of the beam using 
the general expression of Eq. (1), K, and the three models in Eq. (3), Ki, in Eq. (5), 
K2-, and in Eq. (8). A'.-s. 



Ea [keV] 


Ki 


K2 




K 




15 


4.381E+04 


4.045E+04 


7.393E+04 


4.38E+04 


l.OOE+00 


10 


4.073E+03 


3.762E+03 


6.877E+03 


4.11E+03 


l.OlE+00 


5 


1.711E+01 


1.580E+01 


2.892E+01 


1.77E+01 


1.03E+00 


2 


2.615E-04 


2.421E-04 


4.487E-04 


2.85E-04 


1.09E+00 


1.8 


5.038E-05 


7.223E-05 


1.344E-04 


5.62E-05 


1.12E+00 


1.5 


2.339E-06 


3.850E-06 


7.343E-06 


3.34E-06 


1.43E+00 


1.2 


3.613E-08 


7.474E-08 


2.265E-07 


7.84E-07 


2.17E+01 


1 


8.252E-10 


7.711E-10 


5.678E-08 


2.82E-07 


3.42E+02 



2 keV and 1 kcV, we have found that Ki and K2 have still close values, while 
the rate constant K is much larger. It is interesting to note that the constants 
and K have relatively close values. Thus, it is possible to conclude, that, 
for correct estimations of the rate constants, it is quite possible to use the ex- 
pressions shown in (5). The results of these calculations show that the wings 
of the momentum distribution are very important for a correct evaluations of 
the reaction rates. In the last column of table 1 we presented the factor which 
characterizes the deviation of the rate in cause of non-ideal plasmas. 

It is interesting also to estimate the role of the screening effect on the reaction 
rate and to compare it to the considered mechanisms. For that purpose the 
calculations were performed within the framework of the proposed model, but 
with addition of a screening potential Ue = 28 eV. Its influence was taken 
into account in accordance with expression (4). Such value of the potential 
seems realistic for the experimental conditions [1]. The results of calculations 
are presented in the table 2. 

Here we have not presented the results of the more general computations using 
equation (1), because we have obtained rather correct estimates using model 
(5). We find a weak influence of the screening effect, using the reasonable value 
of the potential, which agrees with the results of [1,2]. Taking into account 
quantum corrections, the theoretical evaluations of the rates increase, in the 
low-energy range at 1-2 keV and less, if compared to the rates evaluated 
without the quantum effect. 
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Table 2 

Same as table 1 taking into account the screening effect according to Eq (4) with 
Ue = 28 eV^ 



En [keVl 


± 


Ko 

Z 


^ O 


*■ o / ± 


15 


4.474E+04 


4.133E+04 


7.552E+04 


1.6879937 


10 


4.237E+03 


3.914E+03 


7.152E+03 


1.6880552 


5 


1.911E+01 


1.766E+01 


3.232E+01 


1.69143 


2 


4.059E-04 


5.370E-04 


9.931E-04 


2.4466201 


1.8 


8.433E-05 


1.184E-04 


2.197E-04 


2.6056397 


1.5 


4.605E-06 


7.336E-06 


1.383E-05 


3.0034846 


1.2 


9.315E-08 


1.819E-07 


4.320E-07 


4.638S3085 


1 


2.863E-09 


6.951E-09 


7.386E-08 


25.801192 
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